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Abstract 

We present predictions for exclusive photoproduction of $-meson on proton 
at large transfer, where we use a quark-diquark structure model for the proton. 
Extrapolation from our results to lower transfers is comparable in magnitude with 
available data in that range. This may support the diquark model in its ability 
to provide, for that process, an appropriate link between diffractive physics at low 
transfer and the standard semi-perturbative approach of hard exclusive processes 
at very large transfer, where the proton recovers its three-quark structure. 



I. INTRODUCTION 

Up to now, exclusive photoproduction of a vector meson V {V = p, u, $, J/"^) 
from proton (reaction (1)) have been measured mainly at very low values of t, t 
being the opposite of the squared momentum transfer at the proton vertex. 

7 + p^y + p (1) 

In that region (say, t < 1 GeV^), and in a wide energy range up to HERA 
energies, the observed characteristics of photoproduction of the lighter mesons 
are those of a soft diffractive process. As shown by Donnachie and Landshoff 
im], this can be well described for the most part in a picture using both Vector 
Dominance Model (VDM) and Pomeron phenomenology : there, the incoming 
photon is assumed to convert into a vector meson which afterwards exchanges a 
soft Pomeron with the proton target. In this respect, one may consider reaction (1) 
in this small-transfer range as a good testing bench for Pomeron Physics. Indeed, 
that picture works nicely in case of photoproduction of the lighter mesons p, $ and 
uj 0]. However, it fails to reproduce the energy dependence of the cross section for 
J/\l' photoproduction 0. 

In a QCD-inspired picture, Pomeron exchange is commonly modelled as the 
effective exchange of two non-perturbative gluons. Donnachie and Landshoff im- 
proved substantially that picture 0. When applied to photoproduction of light 
mesons, their two-gluon exchange model leads to results very similar to those pro- 
vided by the Pomeron-exchange model of the same authors . However, it also 
fails to describe the energy dependence of J/^ photoproduction, as well as those 
of virtual photoproductions of p, $ and J/\l', as observed at HERA p. 

Following the works of Ryskin and of Brodsky et al. [0], this may reveal that 
QCD perturbative effects enter the game, since a large momentum scale (either 
the mass of the heavy vector meson produced in case of J/\l/, or the high Q^ of the 



virtual photon in case of p) appears in the reaction. In their approach, Brodsky et 
ah wrote the aniphtude of the process as the product of three terms : an aniphtude 
describing the breaking of the virtual photon into a qq pair, the valence qq wave 
function of the vector meson and an amplitude describing a non-perturbative two- 
gluon structure of the proton. Actually, the latter amplitude plays a crucial role in 
this approach since the dependence on energy of the cross section directly reflects 
the small-x behaviour of the gluon momentum distribution in the proton. In this 
way, one can account for the rapid rise with energy of the cross sections, observed 
at HERA. 

In this paper, we consider (real) photoproduction of $ on proton at larger t. 
Since the $-meson is a pure ss state and the strangeness content of the nucleon 
wave function is probably small, that process is dominated, at lowest order in QCD, 
by the two-gluon exchange mechanism and thus provides a unique way to study 
the latter. This process has been already investigated at moderate t by Laget and 
Mendez-Galain Q, using the non-perturbative picture of Donnachie and Landshoff. 
Their prediction for the differential cross section da/dt at infinite photon energy 
exhibits a node at t = 2.4GeV^, which, in turn, could serve as a test of the model. 
Unfortunately, the only existing data for ^-production with real photons or virtual 
photons correspond to very low values of t. 

Since higher values of t provide larger momentum scales, it is tempting to 
apply in that range the semi-perturbative approach of hard exclusive processes 
developped long ago by Brodsky, Farrar, Lepage (BFL) 0, and by Chernyak and 
Zhitnitsky (CZ) [|T^. This approach has been described and discussed many times 
in the literature, and shown to provide a correct order of magnitude of numerous 
exclusive amplitudes. Indeed, Farrar et al. already applied it to various photopro- 



duction processes JIT . 

We just remind here that in this formalism, the amplitude of a given exclusive 



process is obtained by a convolution formula of relativistic hadron wave functions 
with elementary hard scattering amplitudes involving the valence quarks and an- 
tiquarks of the hadrons taking part in the reaction. When using the baryon wave 
functions derived from QCD sum rules by (CZ), this method provides a correct or- 
der of magnitude for many exclusive amplitudes involving baryons. However, there 
exist important subasymptotic helicity-flip effects that do not fit the above picture 
where any spin effect is to be described by the so-called helicity conservation rule 



12| . While this rule should be valid asymptotically, it appears to be inconsistent 
with most experimental data at intermediate energies. To cure this failure for pro- 
cesses involving baryons, a quark-diquark structure of baryons has been proposed 



13| . In that alternative picture, two of the three quarks of a baryon are cluster- 
ing together in a diquark structure. In the subasymptotic region, diquarks are 
supposed to act as quasi-elementaries constituents having direct couplings with 
photons and gluons. Helicity-flips are then caused by vector diquarks. On the 
other hand, diquarks should asymptotically dissolve into quarks, restoring in this 
way the usual three-quark picture of baryons. The diquark hypothesis provides 
natural explanations for many phenomena that are otherwise difficult to describe 
by standard models fl^ . In the following, we apply that picture in the framework 



of the (BFL) scheme, as a first semi-perturbative calculation of the process under 
study, from moderate to large values of t. 

Another point of theoretical interest in the study of elastic $-photoproduction is 
found in the structure of the amplitude of the underlying hard scattering process. 
Indeed, that amplitude exhibits singularities coming from on-shell quark lines. 
Farrar et al. |T^ have shown that in fact, any exclusive photoproduction process 



is, to leading twist, insensitive to long distance physics and do not require Sudakov 
resummation. The propagator singularities are integrable and their presence does 
not affect the validity of the hard scattering approach. The appearance of an 



imaginary part of the amplitude at leading order in Ug is thus considered as a non 
trivial prediction of perturbative QCD. 

Studies of J/\E' or rjc photoproduction are of course of the same interest as that 
of $ photoproduction, since the charm content of the proton is probably negligible 
too, or even inexistent. However, on one hand, the high value of the c-quark mass is 
a source of computational complications, and, in the other hand, more complicated 
graphs are involved in case of tjc production. So, we leave these two processes for 
future investigations. 

In Section 2 a short description of the quark-diquark model for exclusive pro- 
cesses involving the proton is presented. The details of calculation of the hard 
scattering amplitude for ^-photoproduction is given in Section 3. In Section 4 we 
give our numerical results and concluding remarks. 

II. THE QUARK-DIQUARK MODEL 

Let us first notice that, formally, the amplitudes of the two processes 7 -|-p -^ 
V + p and V -^ p + p + '~^ are just related by crossing. Two of us already studied 
the decay process J/\E' -^ p + p + "^^ using the quark-diquark model structure of 
the proton. So, we largely refer to our previous paper |]l6l for notations. 

The formalism we are using below is the same as that of (BFL), except that the 
three-body structure of the proton is replaced by a two-body one. To lowest order 
in QCD, the photoproduction of $ on proton is thus described by the generic graph 
of Fig. 1. The corresponding amplitude is obtained here too from a convolution 
formula 

[dx] [dx'] [dy] -^T.^^^el^fl^^I^-' (2) 

where [dx] = 6{1 — xi — X2)dxidx2, [dx'] = 6{1 — x[ — X2)dx[dx'2, [dy] = 
5(1 — yi — y2)dyidy2. Here, coUinearity of the constituents with the parent hadron 
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is assumed : Xi = x (resp. x[ = x') is the four-momentum fraction of the quark 
inside the ingoOing (resp. outgoing) proton, and X2 = 1 — x (resp. X2 = 1 — x') that 
of the accompanying diquark ; yi = y (resp. y2 = 1 — y) is the four-momentum 
fraction of the strange quark (resp. that of the strange antiquark) inside the $- 
meson. The tensor T^i^^a is the amphtude for the subprocess gg'~f -^ $ with two 
space-hke gluons having four momenta g = xp — x'p' and G = (1 — x)p — (1 — x')p' 
respectively, p being the four-momentum of the ingoing proton and p' that of the 
outgoing proton ; e?^-) and e? -^ are polarization vectors for the $ and for the photon 
respectively ; I^'^ is a tensor amplitude describing the two-gluon scattering by a 
quark-diquark system ; K is the overall normalization factor : 



Att' 



2. 



K = V47raes f^C (3) 

9v6 

with the color factor C = — 2/(3\/3) and the $ decay constant /^ ~ ISOMeV. 



For the sake of consistency of the model, we have neglected the masses of the 
constituents as well as those of the parent hadrons, whenever possible. This led us 
to use the relativistic form of the $ wave function. Depending on the helicity h of 
the $, it is given by 

^,^ = ^7=^^ss< (4) 

V24 V3,„i„, 1^ ^ ^J) 0^(^) for h = ±l 

where 4>L{y) and (priv) are normalized ^-distributions for, respectively, a lon- 
gitudinally and a transversally polarized meson. 

In that approximation, and due to the particular structure of the amplitude of 
the subprocess (odd number of 7-matrices) , it appears then that only longitudinal 
$ are produced. 



The tensor amplitude T... in (2) is thus simply obtained from the amplitude 
of 37 -^ e+e^ with massless electrons, by removing the coupling constant, using 
appropriate four-momenta and making the substitution 

K+t/e- ^f (5) 

The wave functions of mesons have been derived by (CZ) from QCD sum rules 
technics. We here use their longitudinal- $ wave function 

04i/) = 61/(1-1/) {1/(1 -y)+ 0.8} (6) 



that can be found in [T^, p. 259. 

Allowing for both scalar (S) and vector {V) diquarks, a quark-diquark proton 
state corresponding to a proton helicity "up" or "down" takes on the general form : 



\p^^ > ~ - /s [ 201 (a;) + 03 (x) ] Siud) u^^ ± U [ 02 (x) (v^ V±iud) u^^ 

-2 V±{uu) rf^TJ +03 (x) [V2 Vo{uu) rfTi - Vo{ud) u^^} ] (7) 



where : Vh{qiq2) is an isovector-(pseudo) vector diquark state made of two 
quarks having flavors gi and q2, h = 0, ±1 being its helicity ; S{ud) is the isoscalar- 
scalar diquark state. 

The 0j(x) are normalized wave functions, and /s and /v are normalization 
constants that may be chosen unequal to allow for various admixtures of scalar 
and vector components. Expressions of diquark-gluon couplings have been given 



in refs [|T7|. In the space-like channel {g + D ^ D') these expressions are, using 



obvious notations and omitting color factors as well as coupling constants : 

{S' S), = F, {D, + D'^) (8) 

for a pair of scalar diquarks and 



-Fs(e'};r.D)ie'},.D')iD, + D',) (9) 

for a pair of vector diquarks Q. 

The F above are the diquark form factors depending on Q^ = —G^ = —{D — 
D'Y . A possible parametrization, aiming at describing the natural evolution of 
the diquark model into the usual three-quark picture, has been proposed by the 
authors of ref [17]. It has the following form : 



Fsm = XT^ph^.F,m = x 



-Q^ + QV W + Ql) 

F2{Q') = (1 + K)F[iQ'), F,{Q') = $P^, (10) 



with 



X 



a,mi<^s{Ql) for Q^>Ql 

1 for Q2<Q2 
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fcv being the anomalous magnetic moment of the vector diquark. The value ky,=l 
is commonly assumed. The above-defined evolutionary picture also induces one 
to use, for the sake of consistency, coupling constants of the running form and to 
set Og = as{—g'^)oLs{—G'^)-, with —g'^ = txx', yet restricting as to some maximum 
value Ci. Let us remind that setting the factor x in diquark form factors provides 
the correct power of a^'s in amplitudes at large transfer. 



^We do not consider here a possible mixed coupling involving both scalar and vector 
diquarks, as it is commonly expected to give a small contribution. 



In IjT^, we used such a parametrization to fit the proton magnetic form factor 
Gm in the space-hke region. Modelhng the momentum fraction distributions by a 
wave function of asymptotic form, i.e. taking 

01 (x) = (l)2{x) = Mx) = <Pas{x) = 20x(l - xf (12) 

a quite good fit were obtained with the following values of parameters : 

/s = 40 MeV , /v = 96 MeV , Ql = 2 GeV^ , Ql = 2.3 GeV^ 

for ci = 0.3 (13) 

where, as said above, ci is the maximum allowed value of the running coupling 
constant ctg- Given that success, we used the same parametrization for the present 
calculation. 

III. THE HARD-SCATTERING SUB-AMPLITUDES 

As already mentionned, according to the model here used, longitudinal $ are 
preferentially produced. There are then, a priori, eight dominant helicity ampli- 
tudes describing the process, which we denote by T^^ ^ where A A' and A are the 
helicities of, respectively, the ingoing proton, the outgoing proton and the real in- 
cident photon. Thanks to parity and rotational invariance, that number in fact 
reduces to four and we have : 

(14) 

To be more specific, let us now concentrate on the calculation of the amplitude 
T1^^+. To compute the amplitudes, we have chosen for convenience the polarization 
states of the particles according to a "t-channel helicity-coupling scheme" |jl8|. In 
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that scheme, the photon hehcities represent photon spin projections on the "vertex 
plane" defined by the $ and photon four-momenta, and the photon polarization 
vectors eL-, are then perpendicular to that plane. Similarly, the helicities of the 
protons are projections of their respective spins on the vertex plane defined by 
their two four-momenta. 

From eq. (2) we thus obtain 

T''^ = K^ j[dx][dx'][dy]— f- -0L(y)rTT+ (15) 

t^ J x{\ — x)x'{\ — x') 

where 



TTT+ = 12v/i^F2(Q2)(i _ x){l - x') cot(^/2)x 

1 \ {l-l)A yB y(l-y)Csin\ei2) \ 
y{l -y)\ d' d dd' J 



(16) 



s being the center-of-mass energy squared, 9 the $ emission angle relative 
to the incident photon direction in the center-of-mass frame (sin^(^/2) = t/s), 
Q^ ~ t(l — x){l —x'), fi^ the diquark mass usually taken equal to 600 MeV ; d and 
d' are the s-quark propagator factors : 

d = xx' sin^(6'/2) + y{x cos^{9/2) - x') - ie 

(17) 

d' = {1- x){l - x') sin2(^/2) + (1 - y){{l - x) cos2(^/2) - (1 - x')) - ie 
where, following the usual prescription, e — > 0^. Finally, A, B and C are given 

by 
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A = 03(a:')02(a;)((l - x') sm\9/2) + {1 - y) cos\9/2)) 
-02(x')03(a:)((l - x) sin2(^/2) - (1 - y)) 

B = (l)2{x')(l)s{x){x'sm\e/2)+ycos'^{e/2)) 

(18) 

-03(a;')02(a:)(a:sin2(^/2)-y) 

C = (j)2{x')(j)3{x) {{y{l -x)- x{l - x)) cos2(^/2) + x'{l -y)- x'{l - x')} 
+03 (a;') 02 (a;) {{x{l - y) - x{l - x)) cos'^{9/2) + y{l - x') - x'{l - x'))} 

From eq. (17), it is clear that the kernel eq. (15) has singularities within the 
domain of integration, since the real parts of d and d' have zeroes in x respectively 
located at 

^ x^ ^ . 

^° l/cos2(^/2)+x'sin2(^/2) - 

and (19) 



il-x')il-y) 
[1 - y) cos^{9/2) + (1 - x') sm\9/2) 



^1 = 1- ^ ...... .\. \, . 2^.,.. < 1 



These singularities correspond to one or the two exchanged s-quarks going on- 

shell in the graph of Fig 1. It is important to notice that when x' = y the two 

1 1 

zeroes coincide and are both equal to x' (or y). The one-pole terms (~ ~; or ~ — ) 

can be treated readily, using the general formula 

^ v(-\+m6{u) (20) 



u — le \u 

where V denotes the principal value. The two-pole term ~ — — corresponds 

dd' 

to the graph where the photon line is sandwiched between the two gluon lines. 
Setting r = 3fJ((i),r' = '^{d') (SR means real part), one gets 
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(21) 

It appears that the product of two delta-functions, which is apparently the most 
singular term, leads in fact to a null contribution. This is to be imputed to the 
fact that the amplitude of a fermion-antifermion-vector-meson vertex is zero when 
all particles are massless. 

Let us first consider the imaginary part of the full amplitude. It has the general 
form 

Cr{x,x\y)5{r) + C2{x,x\y)5{r') + C,{x.x\y) \v (i) '^(r') + ^ (^) Kr)] (22) 

that can be trivially integrated by hand over the variable x, yielding an expres- 
sion of the form : 

Ci{zQ,x',y) C2{zi,x',y) 
a a' 

(23) 

l^C,{x,,x',y)V (-) + -Cs{zo,x',y)V (^) 
a' \rJx=zi a \r' J x=zo 

where a = x' sm^{9/2)+y cos'^{9/2) and a' = (l-x') sin^{9/2) + {l-y) cos'^{9/2). 

One must be cautious with the two last terms as they are source of difficulties in 

the subsequent (numerical) integrations, as now explained. Since 



a' \rJx=zi a \r' J x=zo sin^ ^ \{x' - y^ 



3^(r) =r^(3) = -27^ ( 7Z^-T^ 1 (24) 



a double-pole-like term ~ ^/{x' — yY appears. However, that "singularity" is 
tempered by zeroes of the factors C3 when x' = y (C3 oc {x' — y)). These zeroes, 
which are of degree one, are reminiscence of the already mentionned fact that the 
amplitude of a fermion-antifermion-vector-meson vertex is zero when all particles 
are massless ; and, precisely, the two exchanged s-quarks that are coupled to the 
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real photon in the corresponding "singular" Feynman graph get both massless when 
X = x' = y. Consequently, the (seemingly) double-pole reduces to a simple-pole : 

1 \ _ / 1 ~ 



{x' -y)V 



V 



(25) 



{x' — yY j ' \x' — y ^ 

In order to manage this fact in a cautious way, we have proceeded as follows. 
First, we split the products of wave functions into symmetrical and antisymmet- 
rical parts : 



5*23(3;, x') = -{(t)2{x')(t)'i{x) -h 03 (a;') 02 (a;)} 



^23(3;, X') = -{(t)2{x')(t)z{x) - (t)'i{x')(t)2{x)} 



(26) 



Of course, this operation is useful only when 02 7^ 03 and is thus inoperant for 
the symmetrical parametrization used in this paper. We just present it here for 
further applications. The imaginary part of the factor in brackets in formula (16) 
may then be rewritten as 



'k5{x - Zq) {S'23(2o, x')Sq + ^23(^0, x')Aq} + 



(27) 



'k5{x - zi) {S22.{zi,x')Si + ^23(^1, a;')^i} 

In order to maximally reduce the effect of the pseudo-pole oc 1/{x' — y), we 
further made the shift x' = y + {x' — y) and the appropriate simplifications in all 
coefficients S and A. We then arrived at the more manageable expressions : 



So = ^ \x'{x' - y) sin^(^/2) + ya{l + cos2(^/2)) + 
(1 - y) (2a;'(2a;' - 1) + cos2(^/2) [{y - x'){2x' - 1) - 2y^ + y sin^ {9 / 2)]) + 



cos2(^/2) 



{x' + y){2x' - 1) + 2y' 



+ 



(2^ 



yM-yKl-2y)^^±^^^V 



x' — y 
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Ao = sm^{e/2)\ [x'y + sin^{e/2){x' - yf + x' - y + y{l - y) cos'^{9/2) + 



The other factors Si and Ai are obtained respectively from Sq and — .Ao by the 
simple replacement x'-^l — x',y^l— y. 

Let us now turn to the real part of the amplitude. It appears that different 
terms of the same (large) magnitude compensate each other in the domain of 
integration. To cure that new difficulty which causes numerical uncertainties, we 
decided to put all the expression on the same denominator rr' and, again, to 
introduce symmetrical and antisymmetrical combinations of wave functions so that 
the real part of the factor in brackets in (16) takes on the form 

{S^six, x')S' + A23{x, x')A'} V (i) V (^) (30) 



Then, we set u = x' — x, v = y — x', and rewrited the coefficients S' and A' as 
polynomials in u and v. We thus got 

S' = v''Hs + v^H2 + vHi + Ho (31) 



with 
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Hs = -4mcos^(^/2) + 2sin2(^/2)(l + cos2(^/2))(l - 2x') 

H2 = 2u{l - 2x')(l + 2cos^(^/2)) + sin2(^/2)(l + cos2(^/2))(8a;'(l - x') - 1) 

Hi = u^sm\e/2){l - 2x') + 2usin\e /2)x'{l - x')+ (32) 

2u sin2(^/2)(l - Ax\l - x')) + 2u{<ox'{l - x') - 1) 
-2 sin2(^/2)(l + cos2(^/2))a;'(l - x'){l - 2x') 



and 



Hq = ux'{l - x') {(2m - 2x' + 1) sin^(^/2) - 2(1 - 2x')} 

A' = sm^{9/2) {-2v^ + v^{l - 2u - 2x') + 
t;(2x'(l - x') - u^) + ux'{l - x')} 



(33) 



This trick is supposed to moderate both the above-mentionned cancellations 
and the double pole ~ l/(rr'). At this point, let us notice that when f — » 0, one 
has l/(rr') ~ l/{u'^x'{l — x')) whereas the numerators of the amplitudes behave 
like ux'{l — x') so that the double-pole in u turns into a simple-pole (with no end- 
point singularities). Similarly, when -u — > 0, then l/(rr') ~ l/{v^x'{l — x')) while 
the numerators of the amplitudes are ~ vx'{l — x') so that the amplitude is only 
~ 1/v. 

As for end-point singularities due to the denominator x{l —x)x'{l —x')y{l —y) 
coming from gluon and s-quark propagators, they cause no problem since the above 
denominator is cancelled by an equivalent factor contained in the product of hadron 
wave functions. We assume that this cancellation of end-point singularities is suf- 
ficient, namely, that no extra Sudakov form factor that could a priori suppress 
more drastically those singularities is to be implemented. The problem of pos- 
sible Sudakov suppression of end-point singularities is beyond the scope of the 
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present work. It has been discussed in different contexts by various authors to 
whom we refer the reader |[T9|. For our calculations, given the uncertainties in 
the parametrization of the diquark model in its present form, we think it senseless 
to introduce an additional complication with a Sudakov factor, the precise form 
of which for that model and for the process here studied is yet unknown anyway. 
Nevertheless, we here apply a kind of suppression of end-point singularities by cut- 
ting off the dangerous growth of coupling constants asi—g^) and as{—G'^) in the 
end-points region by means of the parameter ci (see above). 

All amplitudes have been treated in the same way. Moreover, to prevent ad- 
ditional instabilities in our numerical evaluations, we modeled principal values by 
the approximate form 

'P(-)-ZT^. (34) 



z/ e^ -|- z'^ 



with e <^ 1. For instance, in real parts we made the substitution 

^ (A) - 2I 2 2'' a (35) 

\rr J e{ + r^ €2 + r''^ 

Our numerical results are presented in the next section. 



IV. RESULTS AND CONCLUSIONS 

We spent a lot of time in searching the best method of integration, and checking 
the stability of our computations. Given the rather simple dependence on y of the 
integrant, we tried to integrate by hand over y first. This is not a good method 
because it induces spurious singularities in the subsequent integration over x and x' 
through denominators like xcos^(6'/2) —x' which does has a zero in the integration 
domain. We also tried a numerical integration in the complex ?/-plane, but this 
also led to intractable instabilities. So, we finally adopted the strategy described 
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in Section 3. We then varied the parameters e down to a value as small as 10~^ and 
even cross-checked our results using two different programs of integration, namely 
a Gaussian quadrature method (RGAUSS) and a Monte-Carlo method (VEGAS). 
It appears that the real parts of amplitudes are the most intractable : the less the 
value of e, the greater should be the number of calls of the integrant. However, a 
value of e between 10~^ and 10~^ seems to provide the best stability. By chance, 
for the particular parametrization here used, the contributions of real parts appear 
to be much less than those of imaginary parts in the whole kinematical range we 
have investigated, by a factor of a few per cent or less. On the other hand, the 
results obtained for the imaginary parts are very stable, being much less sensitive 
to the value of e. 

Thus, the results we are presenting now take account of imaginary parts of 
amplitudes only. 

In Fig. 2, is shown the momentum transfer distribution we obtain in the range 
2 GeV^ - 10 GeV^ and for two values, 5 GeV (CEBAF) and 70 GeV (HERA), 
of the proton-photon invariant mass W. One sees that this distribution is almost 
independent of W. On the other hand, as expected, the distribution exhibits a 
kind of power-law fall-off, like t~^-^ around 3 GeV^ and like t~^-^ around 9 GeV^. 
Actually, we have found that the very simple form 

with A = 94.5 nb GeV^^, B = -0.113 GeV-"^, C = 0.043 GeV-^ provides an 
excellent fit of our results from t = 2 GeV^ up to values of t as large as 15 GeV^. It 
may be useful in a generator program for a simulation of the process. Integrating 
this form between 2 GeV^ and 10 GeV^ yields a cross-section of about 1.5 nb. 

Also shown for comparison in the same figure is the fit of low-t experimental 



data provided by the ZEUS Collaboration pO|. The observed distribution has an 



exponential fall-off ~ exp(-6t) with b ~ 7.3 GeV^^ for < W >= 70 GeV, in 
agreement with the expectation of a diffractive character of the process in that 
range. 

The comparison in Fig. 2 is indeed encouraging for the diquark model since a 
direct extrapolation from our results towards lower values of t nicely compare in 
magnitude with the above-mentioned data. Let us remind here that according to 
the well-known asymptotic constituent counting-rule, one expects a change of the 
observed exponential fall-off of the t-distribution at low t into a power one as t 
is increasing. The diquark model is just supposed to account for that transition. 
One may thus consider the diquark model as a good candidate in describing future 
data at larger t, most probably with more refined wave functions and diquark 
form factors, thus establishing, at least for that kind of process, a link between 
diffractive physics that holds at low t, and the semi-perturbative approach of hard 
exclusive processes one expects to hold at very large t (the proton recovering there 
its three-quark structure). 

Fig. 3 shows s^da/dt as a function of cos(^) where s = W"^. In the pure 
three-quark picture of the proton, that distribution is predicted to be independent 
of s at large s, provided a^, the strong coupling constant, is taken as a constant. 
Obviously, that scaling law does not hold here. Essentially, this is because we are 
using an evolutionary picture of the diquark structure where the a^'s are expressed 
in running forms. It can be easily checked that one recovers the expected scaling 
law, if one takes both the a^'s and the factor x in diquark form factors as constants. 
However, numerical computations show that this works better for W^ > 10 GeV. 
In this respect, W = 5 GeV is not an asymptotic value. 

Such deviations of the diquark model predictions from the asymptotic scaling 
law have also been obtained by Kroll and co-workers in their recent calculation 
on photoproduction of K and K* mesons off proton, using the diquark model too 
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with a similar parametrization pT 



In our numerical calculations, we also tried two other parametrizations of the 
quark-diquark structure of the proton. The first one is that of Kroll et al. in ||22|| . 
We found that the corresponding contribution of imaginary parts of amplitudes 
alone yields a t-distribution that is larger than that of Fig. 2 by more than one 
order of magnitude, which seems to rule out that parametrization since the corre- 
sponding rates look too high, especially at low t. The second parametrization has 
a quark-diquark wave function derived from the asymetrical three-quark proton 
wave function obtained by King and Sachrajda from QCD sum rules p3|. It is 



described in ||T6[. This time, the contribution of imaginary parts of amplitudes is 
much less than that obtained from the asymptotic wave function. On the other 
hand, real parts seem now to contribute much more. But, as said before, the latter 
are difficult to evaluate properly because of instabilities in numerical computations. 
So, we cannot draw any conclusion at the present time about that parametriza- 
tion. At least, this shows that the study of $-photoproduction at intermediate 
t would allow one to discriminate between various models |^. Unfortunately, 



exclusive photoproduction processes are yet largely unexplored at large transfers 
though they possess in that range an undoubted physics potential. We hope very 
much to dispose, in a not too far future, of new data from CEBAF or HERA, at 
least at intermediate transfers. 
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FIGURES 



FIG. 1. A typical diagram for 7 + p^<I>+pin the quark-diquark picture 



FIG. 2. Diquark-model predictions for da/dt vs t ; solid line : W = 5 GeV ; dashed 
hne : W = 70 GeV. Dash-dotted line : fit of low-t data at < IF >= 70 GeV [20]. 



FIG. 3. Diquark-model predictions for s^ da/dt vs cos 6cm ! solid line : IF = 5 GeV ; 
dotted hne : IF = 20 GeV ; dashed hne : IF = 70 GeV. 
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Fig. 1 : A typical diagram for 7+p— ^$+]3in the quark- di quark picture 
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